UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



ILLINOIS COD E M9-l?2 



TITLE 
DURATION 

TARES IBSD 



PROBLEM 



Square Root Factorization 

Input, output, and tape -handling} computing time is 

negligible 

1) Program tape: same for all problems 

2) Parameter tape: 00 9K 

00 pF 00 pF 
00 IF 00 IF 
00 (F+S)F 00 (PfS)F 
00 (P)F 00 (P)F 
26 987N 
where p » number of variables 

P = number of factors previously extracted 
S = number of factors to be extracted on 
this run 

3) Column sums tape: Column sums are punched in the form 
required for D.N. I., scaled down by two decimal 
places; an 'N 1 follows the last column sum} it is 
assumed that each column sum is reduced by unity 
before being punched. 

k) Problem tape: The data is punched in the form required 
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for the D.N. I., scaled down by 10 } an 'N' is punched 

at the end of each column of correlations; the columns 
should be numbered, starting with zero to facilitate 
the locating of columns later. 

This program uses the square root method to find a 

2 
factorization of the given correlation matrix. 



The basis for this assumption is discussed in Wrigley, C. F., and 
McQuitty, L. L. The square root method of factor analysis t a re- examination 
and a shortened procedure. Unpublished manuscript on file at the University 
of Illinois Library, p« 8. 
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Op. Q%t» Ifcis paper gives a full discussion of the method, its advan- 
tages, its disadvantages, and the detailed computational procedure. 
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OWRAITING PROCEDURE 



ccmEsrrs 



1 ) To begin a run tapes are fed to the computer in the 
following order: Program, Parameter, Column sums . 

2 ) The computer then punches ten sexadecimal characters 
in the form OOOABOOCAB 

The two place sexadecimal number AB converted to decimal 
form, say, CD, indicates the column to he "pivoted" on. 
The computer stops after this number is punched. 

3) The column CD is read in. 

M The computer then punches a factor followed by another 
ten sexadecimal digits and stops. 

5) Read in the previous factor (s). The computer stops 
after each factor. 

6) Read in the column of correlations indicated by the 
set of ten sexadecimal characters Just punched. Steps 

k, 5, 6 are repeated until all factors to be calculated on 
the run have been punched. After the last factor the 
computer automatically punches the column sums of the 
residual correlation matrix (in place of another set of 
sexadecimal digits) and then stops on an OF order. To 
begin a second, or later, run, the master tape, a new 
parameter tape, and the just-obtained column-sums tape 
are read in. Operation then proceeds to cycle through 
steps k-, 5, 6 as on first run. Be sure to read in all 
previous factors before each indicated column of 
correlations . 
The capacity of the routine without magnetic storage 
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drum is 250 variables. All factors may be extracted 
if desired* Six factors caa usually be computed in 
a half -hour Machine session*. 

At the end of each session the tapes of the factor* 
that have been obtained during the session should be con- 
solidated onto one tape. This reduces tape-handling on 
the next computer session. In making this consolidation 
one should be careful not to include the pivot numbers 
that have been punched at the end of each factor tape 
on the consolidated tape. 

A pivot may be made on a column other than the one 
indicated by the computer. The steps which are necessary 
to accomplish this are as follows s 

a) Disregard the ten sexadecimal digits which the 
computer has just printed (step k of the regular procedure). 

b) Place a tape in the reader which has the 
following twenty sexadecimal digits 

26 014 00 000 00 (T(FQ) 00 0(PQ) 
where (PQ) is the two-place sexadecimal equivalent of 
the number of the column that you want to pivot on 
next. Remember that the columns are numbered from 
zero . 

c ) Have the computer operator clear R and the 
Order Counter and then read in the tape described in 
(b) with the white switch. 

d) Proceed with step 5 and continue as on regular 
pivot . 

C1ECJLUG ROUTHBB There are several checks built into the above program. 

As long as there are no errors, these codes do not 
affect the operating procedure. At the same time the 
user must be prepared to act if any of the deviations 
from normal operating sequence that are outlined below 
do occur. The deviations, the mistake each indicates, 
and the courses of action to be taken are listed below. 
Program Tape Read- in . A ten-place sexadecimal number 
and some line-feed characters are punched just before 
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the last order on the Program !&pe is read. This output 
indicates an error in the reading of the Program Oape. 

Action: Have the Illiac operator clear the computer 
and read the Program Tape in again. 

Division HMg^. The computer stops on an improper division 
at memory location O^N. Thl. will occur only when the 
diagonal entry is negative in the column of residual 
correlations being used. A residual diagonal entry can 
he negative in case of either a non-gramian correlation 
matrix (e.g., a matrix of tetrachoric correlations) or 
a mistake in the punching of an earlier factor. 

Action; Punch a tape as follows; 
spaces 260K7 00 000 spaces 

Place this tape in the reader, clear ^ and counter, and 
read in on white switch. The residual column sums will be 
punched, ending the computer run . 

SSSHESffl'to residuals . The computer punches (a) the 
residual column sum as Just calculated, (b) 10 line-feed ,.. . 
characters, (c) the residual column sum as reduced from the 
original column sum, and then the computer stops on the 
order « 5 5 N.» This occurs instead of the punching of 
a new factor immediately after the pivot column of 
correlations has been read. One of the following 
mistakes is indicated; 

(1) an error in reading a factor or correlation tape 

(2) not using the indicated column of correlations 

(3) an error in the punching of an earlier factor 
PO a computer error in reduction of column sums 
Action; Repeat the pivot (i.e., repeat steps 5 

and 6) as if the pivot number had just been punched. If 
this check fails a second time on the same pivot, the 
error most likely occurred on an earlier pivot. 

Punch a tape as follows; 
spaces 26 OK? 00 000 spaces 
Place this tape in the reader. Clear B and the counter, 
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and read the tape in on the white switch. The residml 
column sums will be punched and the caaputer will stop 
on an OF order, ending the computer run. 



Coded "by R. J. Twery 

November 195*4- 



Flow Chart for Square Root 
Factor Analysis 



Input orders, parameters, 



and 



ZrW 



1/J 
i 



ij 



+* 



Compare column sums 



r>k 



Test for end of 
run (P+S-s) 



i 



Punch residual 
column suns 
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Punch number of 
column with 
largest sum (Z) 



Read in an 

obtained 

factor 



Increase the p sums 



I 



Test to see if 
all factors have 
"been read in 

1 =T 



4 



Reduce correlations to 
residuals (s) (0) 
r iZ " r 



iZ "& f zj f ij 



Read in the Zth 
column of R 



Check reduction 
to residuals 



X 



I Stops on 20 05N if 
V ~ 1 check doesn't work 



Stretch to factor loadings 

f = r (s) ■*. ( r (s) )" \ 
r is r iZ * v r ZZ ' 



V 



Reduce column sums 

Zr! 6) zjr!!" 1 ' - f. Lf. 
i ij i ij js i is 



Print factor loadings 



>r 



Increase number of factors 
counter (s ) 



-~ «.-H 



Transfer control with the 
order 26 0K7 if residual 
column sums are wanted 
before end of run 



J_! 



.OCATION 



1 



3 
k 

5 
6 

7 
8 



ORDER 

Code X 1 
00 3K 
00 F 
00 6F 
00 6K 
00 250F 
00 250F 
00 22K 

Code N 3 
00 k3K 

Code P 6 

00 177K 

40 ( )F 



50 177F 
26 22F 
26 80F 
00 80K 
hi l4F 
L5 10F 

4o 4f 

iA 8F 
k6 83F 
L5 ( )F 
LO ( )p 
36 86F 

L5 ^F 

J+0 l4F 

1M- 8F 

42 2L 

00 63F 

L5 ^F 

lA 10F 

ho kF 

LO 9F 

36 90F 



from 95F, 100F 



from 153F, I70F 
216F, 218F 

by 99^F 
from 99^F 



NOTES 
Decimal Order Input 



PAGE 1 



M 9 



from 9 



by 993F 5, 12 
by 2 



from 3 



J 



Print constant 
Beginning of storage 

Decimal Number Sequence Input 
Single Column Print 



Read in column sums 



Compare residual column sums 



LOCATION 



9 

11 

12 

13 
14 

15 
16 

17 
18 

19 

20 
21 
22 

23 

24 

25 



ORDER 



L5 4F 
26 8lF 
26 17F 
00 F 
82 12F 
L5 8F 
k2 82F 
L3 12F 
34 94F 
22 96F 

40 ( )F 

50 94F 
26 22F 
22 112F 
32 32L 
kl 5 V 
L* 6"? 
14 14F 
k6 241 
24 19L 
^0 ( )F 

50 19L 
26 22F 

41 4F 
L4 6F 

42 24L 
LO 6F 
L4 13F 
46 25L 
42 25L 
50 ( )F 

TJ ( )r 

L4 ( )F 
40 ( )F 



NOTES 
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from 8 



from 19F 
from 222F 
from 21F 



by 988 
froai 13, 31 



Print number of column with highest 
sum 



(waste) 
from 13 



by 987F 
from 18, 32 



from 59F 



by 18 
by 21 
by 23 
by 23 



Read in the pivot column of correlations 



Read in factors 



Ac emulate sums of cross-products of 
products of factor loadings to be 
subtracted from correlations 



LOCATION 


ORDER 

-* , — _____ 


N °TSS PAGE 3 M Q 


26 


L5 4F 
1 L4 10F 






27 


LO 9F 




. 




32 29L 






28 


L4 9F 








4o 4f 






29 


26 57F 




■ 1 




L5 5F 


from 27 




30 


t li t n-n 

JJT XV/J 








4o 5F 




1 


31 


LO 12F 




. 




34 94f 






32 


24 19L 








4l 4f 


from 15 -~j 


33 


L4 13F 
42 37L 






, * 


LO 13F 
L4 7F 






35 


46 37L 
46 38L 






36 


00 63F 
00 63F 






37 
38 


L5 ( )F 
LO ( )F 
40 ( )F 
L5 4F 


by 35 

by 33 
by 35 


M_ke subtractions, reducing correlations 
to residuals 


39 


L4 10F 
LO 9F 






40 


36 200F 
L4 9F 
4o 4f 




Go to subroutine checking 


41 




reduction to residuals 




26 33L 


_ 




42 


L5 7F 
L4 i4f 


from 214F 




43 


46 44l 









41 15F 







LOCATION 



45 



k6 



vr 



48 



49 



50 



51 



52 



53 



5* 



55 



56 



57 



58 



59 



60 



61 



GRDIR 
L5 ( )? 

50 44l 

2.6 70F 

40 5F 

41 4f 
L4 fF 
46 48L 

42 4% 
L5 ( )F 
66 5F 
S5 F 

40 ( )F 
40 F 
50 F 
7J 16F 
l4 15F 
40 15F 
L5 4f 
L4 10F 

40 4f 
LO 9F 
32 59L 
L4 9F 
22 46L 
00 IF 
00 IF 
L5 l4F 
L4 8f 
46 59L 
00 IF 
L5 ( )F 
92 171F 

41 4F 
L4 7F 
46 64l 
LO 7F 



KOTSS 
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M 9 



by 43 



from 55 



by 47 



by 47 



Find square root of diagonal 
residual 



Stretch residuals to factor loadings 
and form sin of loadings 



by 58 
from 54 



LOCATION 


ORDER 




NOTES p AGE 5 




62 


lA 8f 
42 65L 








63 
6k 


k2 66l 
50 15F 
TJ ( )F 
66 X6F 


by 6l 


Reduce colian sums to those of 
residual Matrix 




65 


4l f 








66 


L5 ( )F 

SO F 


by 62 






61 


ho ( )F 
L5 4F 
IX 10F 


by 63 






68 


LO 9F 
36 190F 








69 


L4 9F 
40 IfF 








70 


22 6lF 


. 








hi 4f 


from I98F 






71 


L4 7F 
k6 72L 


from 76 






72 


L5 ( )F 
50 72L 


by 71 






73 


26 k3F 
L5 hF 




Print factor loadings 




7^ 


lA 10F 
LO 9F 








75 


36 77L 
Lk 9F 








76 


4o 4f 
26 71L 








77 


92 770F 

92 999F 


from 75 


Print N 
Print spaces 




78 


L5 12F 
Lk 10F 




Increase number of factors 
counter 





M 9 



LOCATION 



79 



80 



81 



82 



83 



8k 



85 



86 



87 



88 



89 



90 



91 



92 



93 



9^ 



ORDER 



>0 12F 
LO 11F 
36 87L 
1*1 14F 
iA 13F 
^6 82L 
lH ( )F 
22 59F 
L5 14F 
L4 10F 
14-0 ikF 
LO 9F 
32 86L 
lh 9F 
26 8lL 
26 L 
kl kF 
Lk 8F 
k6 89L 
*l-l F 
L5 ( )F 
50 89L 
26 k3F 
L5 4F 
iA 10F 
LO 9F 
36 94L 
Lk 9F 
^0 kF 
22 87L 
92 770F 
OF F 
00 70K 
Code R 1 
00 20K 



7 

a 



NOiES 
Test for end of run 
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M.9 



from 86 
by 81 
from 6lF 



from 85 
from 80 
from 93 



by 88 



from 92 



from L25F 



Clear locations where sums of 
cross-products are accumulated 



Print residual column sums 



End 



Square Root Routine 



LOCATION 



ORDER 

81 kOF 
kQ ]AF 
22 92F 
00 F 
00 17K 
92 770K 
L5 l^F 
40 F 

82 ^OF 
92 7T0F 
22 91F 
00 57K 
00 63F 
00 63F 
00 63F 
L5 ^F 
26 101F 
00 63F 
00 63F 
00 63F 
26 163F 
00 63F 
00 63F 
00 63F 
L5 ^F 
22 l^OF 
00 190K 

hi 5F 

L5 7F 

iA 5F 

L5 8f 
Lk 5F 
42 5L 
k6 6l 
50 ( )F 



from 90F 



NOTES 



To make a forced pivot 
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from 109F 



Print pivot number 



from 162F 



from 150F 



from ]A8F 
from 9 



Reduce r. a. r 



by 1 



1 LOCATION 


7" 

ORDER 


T~ 


NOTES PAGE 


7 


: 5 


7J ( )F 


by 1 


Adjust for absence of diagonal 




L4 ( )F 


by 3 


entry in column sums 




6 


40 ( )F 


by 4 






i 


L5 5F 








1 ? 


1 L4 10F 
40 5F 


■ 






': 8 


I LO 9F 




. 




• 


32 150F 








9 


22 L 










00 F 
00 16K 


. 


1 




10 


00 F 




1 






00 100 




Scaling constant 




I 


000 








1 


000 
OOOJ 

oo 9871c 




J 







L5 6f 
46 99F 








1 


40 F 
46 94f 








2 


40 7F 
42 57F 


waste 






3 


46 I75F 
L4 9F 


waste 


Form dependent parameters 




V 


40 13F 
L4 9F 








5. 


4o 8f 
46 57F 


waste 






6 


46 58F 
42 82F 


waste 






7 


46 177F 
24 177F 
00 200K 


— 









4i 4f 


from 12 OF 




1 




4l 998F 


-., -- 




1 
i 



M 9 



LOCATION 



ORDER 



NOTES 
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1 


i — 

L4 7F 

46 2L 


1 from 7 


i 

! 


2 


L5 ( )F 


by 1 




1 


L4 998F 






3 


4o 998F 






' 


00 60F 






4 


00 60F 








L5 4f 






5 


L4 10F 




; 




LO 9F 




1 


6 


36 8L 
L4 9F 






7 


4o 4f 








26 1L 






8 


L5 6F 


from 6 






l4 i4f 


! 




..■■?. 


46 12L 




i 


: 


L5 8F 






10 


l4 i4f 








42 12L 






11 


46 17L 






. 


L5.998F 






"12 


LO ( )F 


by 9 




■ 


LO ( )F 


by 10 


. 


■■"""13 


40 IF 

19 13F 






l4 


L2 IF 






. 


36 122F 






15 


L5 998F 
50 15L 






16 


26 43F 








92 189F 






17 


L5 ( )F 


by 11 






50 17L 






18 


26 43F 








20 92F 


,- 





M 9 



Check reduction to residuals 



LOCATION 







ORDER 
00 218K 

:.* 26 if3F 
k\ 179T 

U ( JF I by 1 
L5 175F 
IM- IGF 
40 I75F 
LO 9F 
30 92F 
iA 9F 
26 1L 
00 975K 
Code X 6 

950GN'9*UJEK 
24 999N 



KOTBS 
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I 



Prepare for re -run of pivot in 

ease of error 



- 



Meiiory Sum Cneck 
Sua Check Nuotber 
Control to Do 0. I 



